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ABSTRACT 

We discuss a variety of photometric indices assembled from the uvby Stromgren system. Our aim is to examine the pros and cons of 
the various indices to find the most suitable one(s) to study the properties of multiple populations in globular clusters (GCs) discovered 
by spectroscopy. We explore in particular the capabilities of indices like mi and c, at different metallicities. We define a new index 
64 = (u - v) - (b - y) to separate first and second stellar generations in GCs of any metal abundance, since it keeps the sensitivity 
to multiple stellar populations over all the metallicity range and at the same time minimizes the sensitivity to photometric errors. We 
detecte clear differences in the red giant branches of the GCs examined, like skewness or bi/multi-modality in color distribution. We 
connect the photometric information with the spectroscopic results on O, Na abundances we obtained in our survey of GCs. Finally, 
we compute the effects of different chemical composition on the Stromgren filters and indices using synthetic spectra. 
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1. Introduction 

Thanks to the advent of efficient multi-object spectrographs at 
8m-class telescopes and consequent extensive surveys, in the 
last few years it has become clear that the presence of multi- 
ple stellar populations is ubiquitous in globular clusters (GCs) 
(see, e.g., Marino et al., 2008; CaiTetta et al., 2009a b, 2010a b, 
I2OI Uljoimson & Pilachowskil 120101) . There are at least two gen- 
erations of stars, the second one formed from the ejecta of a 
fraction of the first gener ation (the so-called polluters, see e.g., 
iGratton et al]|200U 1200 4^. This is a crucial piece of information 
on the early phases of evolution of GCs. It see ms to be related 
to their intrinsic, specific formation mechanism (fCarrettal l2006l) 
and may give constrains to theoretical models of formation of 
massive stellar clusters. 

With some delay with respect to spectroscopy, photometry is 
now also indicating the presence of different stellar populations 
(not necessaril y of different age) i n a subsample of clusters (see 
e.g. lAndersonl Il998t iBedin et al.L [2004). We will discuss some 
of its potentialities in the present paper, including strong sen- 
sitivity to abundances of some crucial element (namely N) and 
availability of large samples. However, most of the "hard" evi- 
dence concerning multiple populations in GCs is still based on 
spectroscopic surveys, where the multiple populations manifest 
through a distribution of stars along the Na-O, Mg-Al, and C-N 
anti-correlations, that have been found so far in all GCs where 
abundances of these elements in a large enough number of stars 
have been measured. Typically, modern spectroscopic surveys 
devoted to these studies includ e at best about 100 stars per clus- 
ter (e.g. ICarretta et al.L l2009a l and references therein) with in- 
ternal errors in abundance, derived from intermediate and high 
resolution spectra, typically representing about 8-10% of the to- 
tal variance. While various general information about early evo- 
lution of GCs can be obtained from such distributions (e.g. the 
number ratios between first and second generations and the to- 



tal extent of the Na-O anticorrelation, see Carretta et al. 2009a, 
Carretta et al. 2010c), other require both larger samples and 
smaller internal errors. In particular, this is the case for the dis- 
tinction between smooth continuous distributions and discrete 
ones. Such a distinction is crucial, because it provides a basic 
clue for the nature of the polluters, that m ight be either fast rotat- 
ing massive stars (lifetime a few 10^ vrs: iDecressin et al I I2OO7I) 
or the most massive among intermediate mass stars, undergo- 
ing hot bottom burning duri ng their AGB phase (lifetime a few 
10^ yrs: IVentura et al.Ll200lb . For instance, discrete distributions 
might be only compatible with the AGB polluters, as suggested 
by Renzini (2008). 

Moreover, while multiple populations are an intrinsic prop- 
erty of any GC (so that this migh t be adopted as definition of 
a genuine GC, see ICarretta et al.L [20 10c), the existing statis- 
tics shows also that each GC does show difference in the shape 
and extensions of the Na-O and Mg-Al anticorrelations. In other 
words, apart from a common nucleosynthetic pattern, each clus- 
ter is characterized by its particular history, whose features may 
be studied, in principle, with high resolution spectroscopy of a 
very large number of stars. 

Alternative to spectroscopy, and less time consuming, pho- 
tometry can provide both large samples and smaller internal er- 
rors. The most suitable photometric systems should couple a 
large sensitivity to the abundance variations related to the mul- 
tiple population phenomena, with a good enough efficiency and 
the possibility to be applied to wide-field photometry. Sensitivity 
is obtained by passbands including strong molecular bands of 
CN (bandheads at 3883 and 4216 A) or NH (around 3450 A) or 
CH (around 4300 A), since there are no common use filters sen- 
sible to the (smaller) spectral variations due e.g., to different Na 
abundances. Several photometric systems have been considered, 
including broad band (see e.g . Carretta gt al. 2010d; Lardo et all 
l20Hb . inte rmediate band (s ee lGrundahl et al.Lll999l) . and narrow 
band ones (iLee et al l l2009h . these different methods have advan- 
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Fig. 1. CMDs in yo.Cy o for the nine GCs, seven of which in common with our spectroscopic survey. Blue symbols indicate P stars 
and red symbols IE stars, all confirmed members on the basis of their RVs. For NGC 6205 (=M13) and NGC 6341 (=M92), whose 
spectroscopic data are not homogeneous to the seven other GCs, we use light blue and magenta, respectively; there are two suspect 
AGB stars in NGC 6205 and one confirmed HE st ar in NGC 6341 (note th at the HB is redder than the RGB in these plots). For 
NGC 6397 we display in an inlet the results from iLind etal.l (I2009ll2011h . shown on our same scale in magnitude, and on an 
expanded scale in Cy.o; the separation between Na-poor (light blue) and Na-rich (magenta) stars is perhaps cleaner, but is based on 
a smaller sample (about 30 stars). 



tages and limitations: for instance Sloan filters (IStoughton et al.L 
I200I are very efficient, and high quality photometry can be ob- 
tained even for faint sources. However, sensitivity to the abun- 
dance variations is limited. On the other hand, while narrow 
band photometry may have very high sensitivity, results are not 
well understood yet (see Lee et al. 2009, and Carretta et al. 
2010e). In this paper we will consider in more detail results 
from Stromgren photometry, extending the previous work by 
iGrundahl et all (Il999h and l^ng etall feOOi) . 

Stromgren photometry dStromgreiil Il956l Il966h has been 
originally devised for F-dwarfs and is appUed mostly to the 
study of hot and intermediate spectral type (O to G) main se- 
quence stars. There is however a long series of studies, (see e.g. 
Grundahl et all 1 1 998L 1 1 999l l2002t lAnflionv-Twarog & Twarod. 



20001: lAnth onv-Twarog et al L I2007h~ showing that it can be 
successfully used also for stars on the red giant branch 
(RGB) and on the horizontal branch (HB). In particular, 
lAnthonv-Twarog & Twaroj (l2000i) briefly discussed the effects 
of the "abund ance anomalies " on t he photometry of RGB 
stars in GCs. IGrundahl et al.l (1 1998 1) presented an extended 
Stromgren photometry of NGC6205 (=M13) and noted how the 
RGB stars presented a large spread in the ci index. They at- 
tributed it to the different chemical composition (not different 



iron content) of stars, since the relation between ci and CN, 
CH strength had already been demonstrated (e.g. in M 22 by 
Anthonv-Twarog et al., 1995). A direct measure of N abundance 
was obtained by Yong et al. (2008) and was used to calibrate a 
new index Cy (s ee Sect. 12. 1.11) . A simi lar approach was followed 
for instance by ICarretta et al.l (l2009ah that connected the c, in- 
dex to stars of the first (Na-poor) and of the second (Na-rich) 
generation in NGC 6752. 

ISbordone et al.l (1201 Ih have recently presented a theoretical 
study of the influence of non-standard composition (e.g., en- 
hanced in Na and/or He, and/or light elements like CNO) on 
the stellar models and the resulting isochrones, having in mind 
the particular case of the multipopulations in GCs. They con- 
sidered the cases of Johnson UBVI and Stromgren uvby filters. 
We will discuss later (Sect. 15.1b the similarities and differences 
between their completely theoretical work and our observational 
approach. However, we note immediately that one of the most 
important is that they did all computations for a single metal- 
licit>B value, near the average for the Galactic GC population. 



' Within this paper, the term metallicity is usually adopted as a syn- 
onym of [Fe/H], with an overall abundance pattern characteristics of 
halo stars, that is with an excess of or-elements. 
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Fig. 2. Enlargement of the RGB region for the CMDs in yo,Cyfl of the nine GCs. An arrow shows the assumed magnitude of the 
RGB bump. The Hnes indicate the RGB stars selected to derive average and rm.s. values for Cyfl (see Sect.|4]i. 



while we did instead consider GCs spanning (almost) the entire 
metaUicity range for Milky Way GCs, from about -2.25 to about 
-0.70 dex. 

In this paper we will systematically and deliberately try to 
determine whether photometric data, in particular in Stromgren 
filters, can be used to discriminate between stars of different stel- 
lar generations in GCs, by looking at their RGBs. We see that 
the situation is rather complex, and that there is a strong depen- 
dence on the evolutionary phase and the metallicity considered. 
For instance, the index c,. is efficient in separating different pop- 
ulations in GCs, but only at intermediate/low metallicity, while 
it loses sensitivity for metal-rich clusters (at least for the lower 



RGB). This sensitivity might even change its sign, making the 
interpretation of observational data ambiguous. The mi index is 
instead a good indicator of multiple populations at high metallic- 
ity. Finally, we explore the various indices that is possible to de- 
fine and we select a new index ^4 that uses all the four Stromgren 
filters uvby. We show that it is almost independent from redden- 
ing and that it can be used to efficiently separate different stellar 
generations in GCs over a large range in metal abundance. 

The plan of the paper is as follows: in Sect. 2 we present the 
available data on which we base our analysis and we discuss the 
separation of GC stars in first and second generations as seen by 
the classical indices already known: c,, and mi. In Sect. 3 we ex- 
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Fig. 3. CMD (not corrected for reddening) of oj Cen with stars of 
different [Fe/H] (from Johnson & Pilachowski 2010) indicated 
by different colors: blue stars have [Fe/H]< -1.6, red ones > -1, 
and green ones are in between. 



Table l.List of GCs. 



GC 


Other 


[Fe/H] 


E(S - V) 


yoCRGB/junip) 


Telescope 


NGC 104 


47 Tuc 


-0.76 


0.04 


14.40 


Danish 


JNGC DoJo 


Nil V 


-0.82 


0./5 


14. UU 


JNIJI 


NGC 362 




-1.17 


0.05 


15.30 


Danish 


NGC 1851 




-1.18 


0.02 


16.10 


Danish 


NGC 288 




-1.32 


0.03 


15.35 


Danish 


NGC 6752 




-1.55 


0.04 


13.50 


Danish 


NGC 6205 


M13 


-1.58 


0.02 


14.75 


NOT 


NGC 6397 




-1.99 


0.18 


12.10 


Danish 


NGC 6341 


M92 


-2.25 


0.02 


14.60 


NOT 



Notes. [Fe/H]: Carretta et al. (2009c) and Carretta et al. in preparation 
for NGC 362; E(S - V): Harris (1996); .yo(RGB,„„„p): this work. 



2.1. P, I, and E stars in Stromgren photometry 



amine other possibilities and present our new index ^4, exploring 
its sensitivity to metallicity and to the properties of different stel- 
lar generations. In Sect. 4 we study the correlations between CyX) 
and 64 and cluster parameters. In Sect. 5 we interpret our results 
with the aid of synthetic spectra. A summary is given in Sect. 6. 



2. The data 

We used the publicly available Stromgren photometry for 
nine Galactic GCs collected by Grundahl and cowork- 



ers (iGrun dahl et al.|. [19981 1 19991 l2002h and presented by 
ICalamida et al. (20071). who used them to obtain a new calibra- 
tion of the metallicity index mfl These data were obtained at 
two different telescopes (Nordic Optical Telescope -NOT- on 
Canary Islands, Spain, and Danish telescope on La Silla, Chile), 
as indicated in Table[T] Instruments covering very different fields 
of view (about 4' and 1 1' on a side) were used; sometimes a sin- 
gle field was targeted, sometimes two or more. The observations 
generally comprise the GC center, with the exception of the two 
most metal-rich cluster in the sample , NGC 104 (-47 Tu c) and 
NGC 6838 (=M7 1). We refer to the Calam ida et at] (lOO't) paper 
for references and details on the observations and calibrations. 

We made a selection on photometric errors and sharpness 
(a parameter used to separate well measured stars from galax- 
ies, defects, and cosmic rays), retaining only those stars with 
errors in all the four filters < 0.02 mag and [s harpness] < 0.2 . 
We u sed the 2MASS Point Source Catalog dSkrutskie et all 
l2006h to astrometrize the catalogues through software written 
by P. Montegriffo at the Bologna Observatory. We employed 
the E(B - V) values from Harris (1996) and corrected the 
u, V, b,y magnitudes for reddening using the following relations 
Mo = M - 5.231 X EiB - V), vo = v - 4.552 x E(B - V), 
bo = fc- 4.049 x£(fl-y), vn ^ v-3.277x£(^-- V), with coeffi- 
cients taken from Tab. 6 of'Schlegel etd] ([1998 '). From now on, 
all the quantities, colors, and combination of bands are intended 
as dereddened. For seven of the clusters the reddening is very 
low (E(B - y)=0.02 to 0.05). It is sHghdy higher for NGC 6397, 
but only for NGC 6838 some differential reddening is present. 
Finally, we estimated the magnitude of the RGB bump visually 
inspecting the color-magnitude diagrams (CMDs); no sophisti- 
cated technique was applied, since this information is used only 
to roughly guide our selections (see below). 



^ We downloaded the catalogues from the web page 
[http : //www ■ oa-roma . inaf . it/spress/gclusters . html , 



For seven of the clusters we obtained FLAMES spectra of RGB 
stars in our ongoing survey to study the Na-O anticorrelation 
in a large sample of GCs with different global parameters (see 
Caii-etta et al. 2007a for NGC 6752; Carretta et al. 2009a,b 
for NGC 104 (47 Tuc), NGC 6838 (M 71), NGC 288, and 
NGC 6397; Carretta et al. 2010b, 2011 for NGC 1851; Carretta 
et al., in prep, for NGC 362). We counter-identified stars in 
our spectroscopic samples, containing 50-150 objects each, and 
in the Stromgren catalogues. We generally found about 100 or 
more stars in common, except for NGC 104 and NGC 6838, 
where they are about 15. In the last two cases this is due to 
the fact that our FLAMES fiber configurations are positioned 
around the cluster center, whereas the photometric observations 
are located off the cluster center. We use these stars to try to un- 
derstand whether it is possible to discriminate stars of first and 
second generation (primordial -P- and intermediate plus extreme 
-I plus E- respectively, in our notation: see Carretta et al. 2009a 
for the exact definition) on the basis of their position in color- 
magnitude diagrams (CMDs), experimenting several combina- 
tion of filters. All stars analysed spectroscopically are members 
of the clusters, based on their derived radial velocities and chem- 
ical abundances. 

For NGC 6205 (=M13) and NGC 6341 (=M92) we 
did a similar identification, using literature papers. For 
NGC 6205, a cluster presenting an extended Na-O anticorrela- 
tion (ISneden e t al., 2004; Cohen & Melendez, 2005), we used 
the measures of Na in lJohnson et all (I2005i) for the stars in com- 
mon with our photometry (22 in total, two of which seem AGB, 
not RGB stars) to distinguish P and IE stars, since this paper 
includes more stars than the two cited above. Fo r M92 there are 
no e xtensive studies of the Na-O anticorrelations (lArmoskv et al.l 
1994 measured Na and O in o r ily n ine stars), but 34 stars have Na 
abundances in ISneden et al.l (12000 ). We counter identified ob- 
jects in Sneden et al. with the photometric data and found six 
stars in common (one of them is an HB star). 

2.1 .1 . The c,. index 

Yong et al. (2008) defined the index c,, (see Table |2]for a sum- 
mary of all the usual and new Stromgren indices used in the 
present paper) in their paper on NGC 6752; it traces the N abun- 
dance and was introduced because it is insensible to the tempera- 
ture, at first order, as visible from the near-verticality of the RGB 
at least up to the level of the RGB bump (see Fig. [T]l. The de- 
pendence of Cy to N abundance comes from its definition, since 
it contains both the u filter, where the strongest NH bands are 
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Fig. 4. CMDs in mo, yo for three monometallic GCs and NGC 1851. Blue and red symbols have the same meaning of Fig. [T] darker 
grey points, visible for NGC 6752 and NGC 1851, are blue HB stars (as indicated by their b - y color). 
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magenta for NGC 6205 (=M13) and NGC6341 (=M92)). Note that the HB is redder than the RGB in these plots for intermediate 
and low metallicity and th at the blue HB and the blue stragglers are almost horizontal sequences. As in Fig. 1, for NGC 6397 we 
also show the results from lLind et al.l (2009, 201 1); the separation between Na-poor (light blue) and Na-rich (magenta) stars seems 
cleaner, due to the better data quality and the smaller sample. 



present (at ~ 3450 A) and the v one, where the CN features are 
present (they are weaker but the filter is weighted twice in this 
index). Given the definition of c,., we might then expect sensi- 
tivity to N abundance but also a complex behaviour In fact, at 
low metal abundances and high temperatures only the NH con- 
tribution should be important, CN being negligible due to the 



combination of the quadratic dependence on metal abundances 
(N and C) and the strongest sensitivity to temperature of its for- 
mation. However, at high metallicity, absorption due to CN is 
much more relevant, and may dampen or even change the sign 
of the sensitivity of c,, on N abundances. More details are given 
in Sect. |5] 
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Table 2. Stromgren indices used in the paper. 



Index 


Definition 


Alternative 


Sensitive to: 


m-i 


{v-b)-(b-y) 


v-2xb+y 


metallicity 


Cl 


(m - v) — (v - b) 


u-lxv+b 


gravity 


Cy 


ci-ib-y) 


;< - 2 X V + _v 


gravity & N 


5a 


iu-v)-(b- y) 


C'l + mi 


new index 


S4' 


(64 - 6ab)I{6ar - Sab) 




no T^ff dependence 


6a' 


S'n X (5AR,b - SAB.b) 







Notes. Sab and Sar are the ^4 values derived using the blue and red 
polynomial at the yo of each star (see text and Fig.[6l(. 6AB.b and SAR,b are 
the same, but derived at the magnitude of the bump. 



Fig. [T] shows the 3'o,Cv,o CMDs for the nine GCs, seven of 
which are in the FLAMES spectroscopic survey. As evident also 
from the enlargement of the RGBs presented in Fig. |2] (where 
the lines will be described in Sect.|4|, the structure of the RGBs 
is different in the various clusters: the distribution in c,. can 
be rather uniform and define a more or less tight sequence (as 
in NGC 104 (=47 Tuc) or NGC6341 (=M92)), or present a 
wide dispersion (NGC 1851), or show indications of skewness 
NGC 6205 (=M13) and bi-modaUty (NGC 288, NGC 6752). In 
Figure 1 we identify P and IE stars with different colors (blue for 
P and red for IE stars, and we maintain the color-coding through- 
out the paper). The P stars (i.e., Na and N-poor first generation 
stars) tend to lie on a tight sequence on the blue side of the RGB, 
while the I and E ones (i.e., Na and N-rich stars of the second 
generation), hereafter considered all together, occupy a larger 
color range on the red side. This is most evident in NGC 288, 
NGC 1851, and NGC 6 752, and also clear in NGC 362 and 
NGC 6397. Note that in ICarretta et alJ (|2009a) we could clas- 
sify as P, I, and E only 16 objects for NGC 6397, because we 
conservatively required both Na and O abundances. However, to 
separate in P and IE populations, it is necessary to know only Na, 
so that we have all our 100 stars available for comparison. Only a 
few spectroscopically studied stars are present for NGC 104 and 
NGC 6838, but the two populations do not separate very clearly. 

Su ch a segregation in color has already been found by 
iMarino et al.i (2008) i n M4 (using the Johnson U filter); in 
NGC 6752 by Carretta et all |2009a!), where we identified the 
N-rich stars in Yong et al. (2008) with the Na-rich ones; and in 
NGC 6397 bv lLind et all (1201 ID . using the Stromgren index c,. 
(and tightening the connection between abundance variations of 
different element s). The same effec t (and the same explanation) 
was explored bv lLardo et alJ (|201 ih in nine GCs using the Sloan 
u filter. 

Our first conclusion is therefore that Cy^ seems to efficiently 
separate sequences of first and second generation stars at inter- 
mediate to low metallicities (indicatively for [Fe/H] ^ - 1.3), 
but does not work well for metal-rich clusters like NGC 104 or 
NGC 6838, or even NGC 362. We will see in Sect. |5] that the 
effects of chemical abundances are complex. 

In NGC 1851, even if the two stellar populations, P and IE, 
are segregated in Cy o along the RGB for yo > 15.5, the branch 
has a larger spread in Fig. [T] with respect to the other GCs. 
Recently Carretta et al. (2010b) presented the hypothesis that 
NGC 1 85 1 is the result of a merger of two originally distinct GCs 
with slightly different metallicity, each one with a Na-O anticor- 
relation. The appearance of the plot would be explained if we are 
seeing the effect of a (small) spread in [Fe/H] and likely overall 
CNO content (Yong & Grundahl, 2008), with Na-0 (hence N) 
differences in both the metal-rich and metal-poor components. 
That Cy measures both N and the overall metallicity is imme- 



diately visible using the template cluster with a dispersion in 
[Fe/H]: lj Cen. In Fig.|3]the most metal-poor population defines 
the right edge of the RGB, with stars of increasing metallicity 
populating sequences at more and more negative values of Cy 
(we used the stars studied by Johnson & Pilachowski 2010 with 
Stromgren photometry by Calamida et al. 2007). The spectral 
synthesis calculations reported in Sect.|5]show that a variation of 
[Fe/H] from -2.23 to -0.63, roughly the range observed in ai Cen, 
causes changes in Cy of 0.1 1 mag for subgiants and 0. 18 mag for 
stars at the RGB bump, c,. being smaller in more metal-rich stars. 
Such differences in Cy are comparable to those expected between 
N-rich and N-poor stars. Similar differences are also produced 
by variations in the CNO content, within the limits that might be 
expected for stars in NGC 1851 and cj Cen (I Yong & Grundahli 
I2OO8I; IMarino et al.L I2OI lb . Disentangling these effects can be 
difficult without further information. 

2.2. The mi index 

The Stromgren index usually employed to trace the metallic- 
ity is nil. In Fig- |4]we display the effects of metallicity using 
the mi index, or better its dereddened version, mo, displaying 
three representative monometallic clusters and NGC 1851. In 
NGC 104 (=47 Tuc) the RGB has a spread in mo, but this clus- 
ter is well known to h ave a negligible dispersion in [Fe/H] (e.g., 
ICarretta et aI.L [20Q9d). Thus, this diagram indicates that the in- 
dex mo is sensible also to the abundance of elements other than 
iron at high metallicity (e.g., N is present through the CN bands 
visible in the v filter, see discussion on c, above and Sect. |5]l. 
This dependence is almost not detectable for NGC 362, similar 
in metallicity to NGC 1851, and is undetectable in the interme- 
diate metallicity NGC 6752 (where the apparent spread of the 
RGB fainter than 3'o~13.5 is due to the blue HB stars, which 
form an almost vertical sequence). 

The P and IE stars occupy the bluer and redder parts of 
the RGB, respectively for NGC 104 (and NGC 6838=M71, not 
shown here), while the separation is much less clear in NGC 362 
and is not seen at all for NGC 6752 (and the other clusters of sim- 
ilar or lower metallicity). In NGC 1851 there is also a spread, not 
imputable to the blue HB or simply to a metallicity difference, 
since the dispersion in [Fe/H] determined from spectroscopy 
should yield only a spread of 0.01 mag in the nii index, accord- 
ing to our spectral synthesis (see Sect|5]l. This cluster is very 
comp lex and many of its feature are discussed in (ICarretta et al.L 
I2OI ll) . Our second conclusion is then that mo is a decent indi- 
cator of first and second generation, but only at high metallicity 
(indicatively for [Fe/H]^ — 1.0). 

To summarize the results of this discussion, the index c,. is 
a good indicator of different N abundances (i.e., distinct stellar 
generations) in metal-poor clusters. For metal-rich clusters, like 
NGC 104, the variations due to N in the v-b and u-v colors are 
similar and the net effect is to cancel or even to change the sign 
of the sensitivity of Cy to different populations; in this metallicity 
regime nii is a better indicator 

3. A new index, 64 

To reduce as much as possible any degeneracy with metallicity, 
we test here a new Stromgren index, that we call 64, defined as 
64 = (u - v) - (b - y), or 64 - ci + mi. We think that it can be 
an useful indicator for all metallicities, and that it minimizes the 
dependence on photometric errors, since it does not weight any 
of the photometric filters twice. Given the limit to photometric 
errors to < 0.02 mag adopted to select stars, the errors in 64 
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Fig. 6. CMDs in 64,yo for the RGBs of the nine GCs, with the RGB bump level indicated by an arrow. The blue and red curves, 
stopped one magnitude brighter than the RGB bump, where the statitics begins to be too scarce, are polynomials describing the 
populations, drawn according to the distributions of P and IE stars in Fig.|5] For NGC 6752 the distribution of stars is trimodal (see 
text and next figure) and an intermediate (green) line is drawn. 



should be < 0.08 mag. However, a more realistic estimate is 
an r.m.s of ~ 0.02 mag. Furthermore, both nii and especially 
Cy do have a stronger weight on the bluest filters, the ones most 
probably having the largest photometric errors, especially for red 
stars. In our new index the photometric errors in (u-v) and (b-y) 
are not correlated with each other. 

We have used in the present paper only dereddened val- 
ues; however 64 would not depend much on reddening any- 



way. Taking from IStromgrenI d 1 9661) the dependence on redden- 
ing for ci and mi, we have that £(64) = 0.02 x E{b - y) = 
0.014 X E(B - V), i.e., ^4 is approximatively reddening-free, at 
first-order approximation. 

To test the ability of ^4 to separate stars of different genera- 
tions, we plot in Fig. |5] the 64,yo CMDs of the nine GCs, using 
the same symbols of Fig.[T] the figure shows that there is a good 
separation between P and IE stars. Also using this new index we 
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Fig. 7. Distribution in ^4" for the nine clusters, showing the large difference among the GCs. The stars are colored in blue or red 
(and green, only for NGC 6752) according to the separation made in 64 (see previous figure); this division is only indicative. The 
colored histograms refer to the blue/red populations. The black, heavy lines represent the cumulative functions of all stars below the 
RGB bump, indicated by an arrow. These cumulative functions (to which the left y-axes refer) are the ones to consider to determine 
whether/where the distributions change slope because of a separation in color (i.e. population). The most notable cases are NGC 288, 
NGC 6752 (and maybe NGC 104 (=47 Tuc) and NGC 362). 



see some structure in the RGBs, most evident for two clusters, 
NGC 288 and NGC 6752; there the distribution of stars in 64 
seems bi-(multi)modal and not continuous. This is best shown 
by Fig. |6] where we display only the RGBs of all nine clus- 
ters on an enlarged scale. NGC 288 appears to have two dis- 
tinct branches, one populated by P, the other by IE stars (see the 
previous figure), while NGC 6752 displays a tri-modal distribu- 
tion (see also Fig. |7]). If we trust the link between 64 and stel- 



lar generation, we could conclude that in NGC 288 there have 
been two episodes of star formation from gas well homogenized 
in each phase, without any appreciable continuous dilution from 
pristine gas, as invoked by models (e.g.. lPrantzos & Charbonnel 
2006) to explain the run of the observed Na-O anticorrelations. 
In NGC 6752 there are three populations, even if their separation 
is not so strikingly evident as in NGC 288. 
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Fig. 8. Relation between ^4" and the Na abundance for the seven GCs in our spectroscopic survey (and NGC 6205=M13, shown 
with different colors, like in Fig. [TJ. Blue and red symbols are for P and IE stars, respectively (the definition is based only on Na 
abundance, to increase the sample sizes). 



The variation of 5^ seems to depend on magnitude, i.e. on 
temperature (Figs.|5]|6]l and also on metallicity, see next Section. 
Ignoring this effect would produce a blurring of the 6\ values and 
there would be the risk of potentially hiding some real features. 
Guided also by the distribution of P and IE stars that we see in 
Fig. |5] we traced reference sequences in the 54,3'o plane. First, 
we fitted a cubic polynomial in yo to the RGB ridge line; second, 
we labeled as "blue" the stars with negative residuals and "red" 
those with positive residuals with respect to it (P and IE stars, 
respectively); and finally we drew the best fitting cubic poly- 
nomial in yo to the two distributions of blue and red stars. The 
results are shown in Fig.|6](where NGC 6752 is actually divided 
in three distributions). 

These curves can help us to eliminate the curvature of the 
RGB using both the blue and the red sequences, and immedi- 
ately compare the distributions of red and blue populations. We 
then define a new index, called d'^ to rectify and normalize this 
index. Of course, with this procedure we risk to normalize to the 
errors in those cases where the intrinsic spread is small with re- 
spect to observational errors (like e.g. the case of NGC 6397). 
We then defined a second index (5^' which is 6'^^ multiplied by the 
r.m.s. scatter around the best fit line for the whole sample (see 
Table 12] for definitions). 6'^ takes into account the fact that the 
blue and red sequences can be more or less separated in color 
in different clusters. This index can be used to better identify 
the presence of multiple distributions along the RGBs. Note that 
given its definition, photometric errors on 6'!^ should be similar 
to those in ^4, with typical rm.s values of ~ 0.02 mag. 



We also note that 5\ and 5'^ are only defined within a lim- 
ited range in magnitude for each cluster, and extrapolation of 
the best fit lines out of these regimes produce results that are 
clearly wrong. Our following discussion only considers those 
stars fainter than the RGB bump (or slightly brigther). 

3. 7. 6\ and the separation of star generations 

Fig.|7]shows the RGBs of the nine GCs "rectified" by using this 
new index, and allows to appreciate that strong differences do 
exist among the GCs in our sample. This division is rather arti- 
ficial in most cases: we simply divided second generation (red- 
der) and first generation (bluer) stars by choosing the stars that 
lie near the red and blue lines in Fig. |6l The clear exceptions 
are NGC 288 and NGC 6752, since the distribution of giant 
stars is clearly bimodal in the first one, and tri-modal in the 
latter NGC 104 (=47 Tuc) and NGC 362 show smaller indi- 
cation of separation, although a bi-modal distribution is not ex- 
cluded; NGC 6205 (=M13) has perhaps the larger fraction of 
IE stars; and NGC 1851 has the largest dispersion in 6'^, possi- 
bly in agreement with the hypothesis of an origin by merger of 
two clusters with different chemical composition (Carretta et al. 
2010b, 2011). 

To effectively see whether we can really discern a separation 
of populations, in the same figure we also present the cumulative 
distribution in 64' (limited to stars below the RGB bump). A 
change in the slope is an indication of separations into discrete 
distributions; NGC 104, NGC 288, and NGC 6752 are the only 
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three GCs where such a feature is clearly present (and maybe 
NGC 362). 



Only in the most evident cases of NGC 288 and NGC 6752 
we estimated the fractions of second (red, red plus green in the 
latter) and first generation (blue) stars as defined by the photom- 
etry and the 64 (and 64") index and compared them with the cor- 
responding fractions of P and IE stars, defined spectroscopically. 
In Carretta et al. (2009a) we found that for all GCs the fraction 
of P stars is more or less constant, at a level of about one-third. 
Here we find the same result for NGC 6752 (blue stars=27%, 
red stars=73%, compared to P=27% and IE=73%). However, 
for NGC 288 the fractions of P and IE stars are reversed (62% 
and 38% from the photometry compared to 33% and 67% from 
the spectroscopy). The split in the RGB of NGC 288 is very 
clear (and is present, in the same sense, also in the Cy_o,yo CMD). 
There are several considerations to be taken into account. First, 
the discrepancy is reduced if we classify stars in P and IE us- 
ing only Na abundances, since in this case we have fractions 
of 42% and 58%, respectively. Second, in NGC 288 the spec- 
troscopic sample is limited, perhaps more than in other GCs, 
to stars significantly brighter than the bulk of the photometric 
sample, while the stronger concentration of blue/P stars is at the 
fainter limit of the RGB, where the actual separation is less ev- 
ident. Third, we see in Figs.[T]and|5]that, at the fainter limits of 
the spectroscopic sample, there are a few IE stars that overlap in 
color the P ones (the reverse does not happen). In this particu- 
lar cluster the division in color/index is quite natural, while the 
spectroscopic classification i n P and IE res t s on th e less defined 
minimum value for Na (see ICarretta et al.L l2009ai for a defini- 
tion); indeed, a slightly different separation in Na would bring 
the fractions of P and IE stars more in agreement with the pho- 
tometric ones. In conclusion, the result we found on the relative 
importance of first and second generation stars (one third and 
two thirds) is statistically valid on the whole sample of GCs, but 
in some clusters there may be individual variations. While this 
issue merits attention and further considerations, this is besides 
the scope of the present paper. 



Since the separation between populations mostly rests, spec- 
troscopically, on the Na abundance, we plot in Fig. [8] the rela- 
tion between 64" and the Na abundance, which is rather tight. 
Here we have divided stars in P and IE using only Na abun- 
dance, to increase the samples (we checked that this approach 
does not generally make a significant difference with respect to 
the criteria for PIE classification as adopted in Carretta et al. 
2009a) and we also restrict to stars below the RGB bump since 
the separation between the sequences, hence the definition of 
64", is better We also show the case for NGC 6205 where, how- 
ever, most of the observed stars are very bright. Nothing can be 
said for NGC 104, and NGC 6838 (=M71) since we have too 
few stars. For the other GCs, there is a very good correlation 
between the 64" index and Na, especially for the intermediate 
metallicity ones (the sensitivity is worse at very low metallicity, 
like for NGC 6397, where abundances are derived from rather 
weak lines). The relations could be interpreted as continuous, 
but we see the possibility of separating stars in two groups, again 
most clearly in NGC 288 and NGC 6752, but also in NGC 1851 
and especially in NGC 362 (where neither 64" or Na show well 
defined partitions). It is sounder, and more natural, to separate 
the populations using both photometric and spectroscopic infor- 
mation, whenever possible. 
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Fig. 9. Relation between metallicity and the value of 64 at the 
reference magnitude My = (i.e., at about the HB level) com- 
puted for the whole RGB and for the blue and red stars (P and 
IE, respectively). 



3.2. 64 and metallicity 

Both nil and c,, have a (complex) dependence on metallicity 
([Fe/H] and N abundance). To explore what is the dependence 
of ^4, we have measured the value of ^4 at a reference magni- 
tude (Mv ^ My - 0, i.e., near the level of the HB) and for the 
different populations, computing three 64 values (blue, average, 
and red), guided by the natural separations and the lines drawn 
in Fig.|6] Given the rather good correlation existing between 64" 
index and Na abundances, we label the blue and red lines as P 
and IE respectively, even though it should be clear that the defi- 
nition of PIE groups should be based on spectroscopy rather than 
on photometric data. 

The values we obtain by this procedure for the clusters we 
examined are plotted in Fig. |9] We see that 64 has a strong 
dependence on metallicity. The whole range of 64 observed in 
NGC 6752 can be obtained if the metallicity varies by about 0.3 
dex (i.e., the difference between NGC 6752 and NGC 288). We 
then expect some difficulties in separating N-rich and N-poor 
populations in clusters like NGC 1851, where there is an inter- 
nal spread in metallicity (Yong & Grundahl 2008 and Carretta 
et al. 2010b find a scatter of about 0.08 and 0.07 dex, respec- 
tively, well above the expected uncertainties). We conclude that 
64 is a very good estimator of the N abundance if a cluster 
is monometallic, as almost all are, at least at the precision we 
can reach with the present day spectra and analysis methods (in 
Carretta et al. 2009c we showed, using UVES spectra of 19 GCs, 
that the internal homogeneity is at about 10%, given an upper 
limit to the scatter of iron of less than 0.05 dex). Where this is 
not valid, as in NGC 1851, it is better to separate the depen- 
dence from [Fe/H] and N using, e.g., mi. Of course, 64 also has 
a dependence on temperature, as seen from the curvature of the 
RGBs. 



4. Cj 0, ^4 and cluster parameters 

We have explored the possible relations of the Cy o and 64 in- 
dices with global cluster parameters. For Cy 0, we concentrated 
on the stars below the RGB bump where the distribution in Cyfl 
is almost vertical, i.e. independent of the temperature (Fig. |2]l. 
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As already discussed, the distribution in Cy o is different from 
cluster to cluster and may present evidence of structures. These 
findings, given the connection between this index and N, must 
be indicative of different distribution in N-richness among the 
GCs. 

For each GC we chose only stars in a two-magnitude bin be- 
low the RGB bump (i.e., y^ump -1-0.25 to ytn„„p+ 2. 25) and selected 
as RGB those falling within 3cr from the average Cy.o value (af- 
ter a K-sigma clipping); the selected region is indicated in the 
figure. The adopted criteria allow us to deal with statistically 
significant samples of stars and the verticality of the RGB in this 
plane allows us to separate rather easily the possible contamina- 
tion by field stars. We measured the average value of Cy.o in the 
defined interval and its r.m.s. scatter. In principle, we expect that 
a higher average value does indicate a higher relative weight of 
the second generation, while a larger r.m.s. scatter could be the 
signature of a larger contribution of the nuclear processing by 
the first generation polluters to the yields budget in the cluster 

We show the best correlations we found for the average Cyfl 
and its r.m.s. in Fig. [TO] As expected and also already visible 
from Fig.|2] the average depends on metallicity (left panel in sec- 
ond row of Fig. [Toll. We have tried to see if there are other signif- 
icant correlations with several structural or chemical properties, 
but results are not very robust. The average increases somewhat 
with the cluster mass (represented by its proxy, the total absolute 
visual magnitude My), and it is larger where the Na-0 anticor- 
relation is more ex tended (as indicated b y the interquartile ra- 
tio IQR([Na/0], see lCarretta et alj|2009d) . Similarly, it is related 
to 5Y, the am ount of He-enhancem ent that can explain the long 
blue HB tails dGratton et al.Ll20l" ^ and which is connected with 
the extension of the Na-O anticorrelation. Some correlation exist 
with the relaxation time at half-mass radius, once 47 Tuc (which 
has the longest relaxation time in the sample) is excluded (right- 
most, upper panel of Fig.fTOb. A full cluster relaxation occurs in 
at least 2 -3 relaxation times and considered the age of 47 Tuc 
(see, e.g.. lGratton et al.L ,2010) this GC should be fully relaxed. 

We note that 47 Tuc was alre ady indicated bv lCarretta et alJ 
( 120 lOd) and iGratton et all ( 1201 Oh to be a pecuhar cluster, with 
rather small values of IQR[Na/0] and He spread with respect 
to other clusters of similar total mass. The same is apparently 
visible in the correlati on of the average C y.o index and My in the 
present data. Recentlv. lLane et al] (1201 Oh suggested that also this 
cluster might be the result of a past merger of two originally dis- 
tinct clusters. However, we must point out that for 47 Tuc the 
number of stars is smaller than in the other GCs, given the more 
external region sampled by the available Stromgren photometry. 
It is possible that at such larger distances from the cluster cen- 
ter, the ratio first-to-second generation stars is shifted toward the 
primordial pop ulation, as expected on theoretical grounds (e.g., 
2008) and a s seen e.g., from the distribution of 
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CN-rich stars ( Norris & Freemanll 19791) . interpreted as Na-rich, 
second-generation ones in the present-day framework. 

The last row of panels in Fig. [TOl shows the best correlations 
we found using the dispersion in the Cy.o index. In the case of the 
r.m.s. scatter, the c orrelations with the spread of He SY (from 
iGrattonet £0120101) and IQR([Na/0]) are significant at more 
than 99%; the first is muTored by the anticorrelation between 
the spread in Cy pa nd the Mmm, the min imum mass reached along 
the HB (also from lGratton et al.ll2010h . 

These findings represent another indication that Cy is tied to 
the different generations through its sensibility to N (and the link 
with Na, O, and He abundances). 

A similar approach can be taken for ^4, and we show in 
Fig. [TT] the results. Also in this case we used stars below the 



RGB bump and determined the average and rm.s. (we used here 
the median, given the less uniform distribution in this index) of 
^4" for the nine clusters. The indications are similar to the previ- 
ous case for My, IQR, 5Y, and M„„„, while the correlation with 
the relaxation time disappears. 



5. Comparison of observed colors with synthetic 
spectra: interpretation 

In order to better understand the observed runs of the classical 
and new indices we obtained from the Stromgren system, we 
computed a number of synthetic spectra over the wavelength re- 
gion of interest (approximately, A - 3000 - 6000 A) for some 
cases of interest. The synthetic spectra were computed using the 
iKurucd ([T99I set of model atmospheres (wiflijhe overshoot- 
ing option switched off), and line lists from lKurucd (Il993h CD- 
ROM's. It is well known that this line lists, while very extensive 
and useful, still have important limitations. In addition, both the 
model atmospheres and our syntheses include approximations 
(e.g. mono-dimensional atmospheres, scattering treated as an ad- 
ditional absorption, etc.). For this reason we only consider dif- 
ferential effects between stars that should be very similar, except 
for details of the chemical composition. 

All synthetic spectra were computed with a wavelength step 
of 0.02 A. Then they were convolved with a Gaussian with a 
FWHM of 2 A, and rebinned at about 1 A sampling for an easier 
handling and display. The sp ectra were fi nally convolved with 
the transmission of the Mvfey dStromgreni 1 1956) and Ca filters 
dAnthonv-Twarog et al.Lll99lh. co llected in the Asiago Database 
on Photometric Systems dMoro & Mun ari. 2000). 

For all choices of atmospheric parameters, we considered as 
reference the spectra computed with the typical abundance pat- 
tern observed in field metal-poor stars, that is [C/Fe] = [N/Fe]=0, 
[O/Fe]=[Mg/Fe] = [Ca/Fe]=0.4, and [Ti/Fe]=0.2. Spectra com- 
puted with these parameters are labeled "N-poor ". They should 
mimic the composition of P stars according to ICarretta et al.l 
(12009a'). For the same choices of atmospheric parameters, 
we computed spectra for atmospheres having [C/Fe]=-0.2, 
[N/Fe]=-i-1.3, and [O/Fe]=-0.1 (and the same abundances for 
Mg, Ca, and Ti). We labeled this second group of spectra "N- 
rich"; they should mimic the spectra of I stars|_| In addition, we 
also explored in a more limited way the influence of C increase, 
computing spectra for [C/Fe]=0.2, [N/Fe]=0.0, [O/Fe]=0.4 ("C- 
rich,N-poor" spectrum) and [C/Fe]=0.2, [N/Fe] = 1.3, [0/Fe]=- 
0.1 ("C-rich,N-rich" spectrum). All these compositions are also 
indicated in Table [3] 

First we considered pairs of "N-poor" and "N-rich" spec- 
tra for the four metallicities [Fe/H]=-2.23, -1.63, -1.23, and 
-0.63. They have effective temperature Teff and surface grav- 
ity log g appropriate for stars at the base of the subgiant branch 
(SGB), at the RGB bump, and at My = -1 (as read from the 
isochrones of log(age)=10.10 bv iBertelU et"al]l2008h . Table [3] 



^ We do not expect any significant change between I and E stars in 
this respect, insofar the abundance of O is larger than that of C. The 
main difference between I and E stars is in the O abundance, which 
changes from [O/Fe] about -0. 1 to about - 1 , but even in this case spec- 
troscopic observations do not indicate the presence of very strong C- 
bands, which are expected as soon as the abundance of C exceeds that 
of O. The filters we are considering do not contain OH bands (e.g., the 
strong ones near 2700 A) but only N molecules. However, N is already 
about 90% of the whole sum of CNO elements in I stars, and the small 
further increase expected between I and E stars (about 5%) does not 
influence our computations. 



11 



E. Carretta et al.: Stromgren photometry and MPs in GCs 



r =-0.40 
r =-0.50 H 



-0.1 



-0.2 - 



-0.1 



-0.2 - 



= 0.43 - 
= 0.53 - 



-0.1 



-0.2 



= 0.22 
= 0.10 



-0.3 



-0.3 



Mv 



0.5 1 
IQR[Na/0] 



-0.3 



2 3 
, (Gyr) 



r =-0.50 
r =-0.58 



-0.1 



-0.2 - 



-0.1 



-0.2 -• 



r =0.60 
r =0.48 



-0.3 ^ 

0.05 pr 

0.04 r 

°^ 0.03 r 

g 0.02 r 

0.01 r 

L 





-2 -1 
[Fe/H] 



r^ = 0.6B . 

I I I I I I I 
0.5 1 
iqR[Na/0] 



-0.3 ^ 
( 

0.05 rr 

0.04 r 

J"- 0.03 r 

g 0.02 r 

^ - 

0.01 r 




0.05 
5Y 



r,= -0.77 

r =-0.71 



0.1 



0.4 0,5 0.6 0.7 0.1 
M„,„ 



0.05 
0.04 
0.03 
g 0.02 
0.01 




r^ = 0.B7 
r =0.82 



0.05 

6Y 



0.1 



Fig. 10. Correlations of the average and r.m.s. for Cy o with some interesting parameters; the two metal-rich clusters are indicated by 
red points, rs and rp are the Spearman's rank and Pearson's linear correlation coefficients, respectively. The only highly significant 
correlations (significance of the Pearson coefficient larger than 99%) are with metallicity for the average value of Cy o and with 
lQR([Na/0]), M,ni„ and 6Y for its r.m.s.. 



lists the values of the atmospheric parameters we used and the 
offset in magnitude between the "N-poor" and "N-rich" spectra 
for the uvby and Ca bands and for various photometric indices 
(mi, ci, Cy, and 64). 

Figure [T2]presents an example for the synthetic spectra ap- 
propriate to stars at the RGB bump (panel a) and above the RGB 
bump (panel b) for [Fe/H]=-1.23; the transmission of the fil- 
ters is also shown. This figure visually demonstrates what has 
already been discussed in the previous sections. The differences 
between the "N-poor" and "N-rich" spectra are essentially due 
to the different strength of the molecular bands: the NH band at 
~ 3400 A, which is stronger in the "N-rich" spectra and falls 
within the u band; and the CN violet band at ~ 4216 A, which 
falls within the v band. The CH A-X band at ~ 4320 A is also 
different between the two spectra, but is near the edge of the v 
band, so that its impact on the Stromgren colors is very small. As 
a consequence, the flux predicted for both the v and the u bands 
is smaller for the "N-rich" spectra than for the reference "N- 
poor" ones. The difference is larger for the u band, where it may 
be as much as 0.2 mag, save for the coolest and more metal-rich 
stars, where the effect is larger in the v band. This is due to both 



saturation and different temperature sensitivity of the molecular 
bands. A sketch of the variations in the u,v,b,y,Ca filters and 
in the three indices mi, Cy, and 64 is presented in Fig. [T3] where 
we can appreciate their different sensitivity to CNO abundances 
along the various evolutionary phases and their run with metal- 
licity. Table[3]and Fig.[T3]can be used to interpret what we have 
seen in the observed CMDs in term of different chemical com- 
position (in other words, of different stellar population). 

While the largest variations between first and second gener- 
ation stars are expected for N, also the C variations can have an 
effect on the photometry. We computed only a couple of spectra, 
and Fig. [12] shows also the results a for C-enhanced composi- 
tion. The corresponding differences in the various filters/indices 
are indicated in Table [3] From these computations we see that, 
for N-poor stars the effect of an increased C abundance is small, 
since the CN bands are not very strong and O is much more 
abundant than C in any case. On the contrary, for N-rich stars, 
an enhancement in C produces large changes to the Stromgren 
magnitudes, especially for the u and v filters. This means large 
differences in indices like c,, and mi (while the differences al- 
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Fig. 11. Correlations of the median value and r.m.s. for ^4 with some interesting parameters. Symbols are as in previous figure. 



most cancel out for 64, since it is the sum of the previous indices, 
whose variations have opposite sign). 

The combination of N and C variations could then explain 
some peculiar cases, like NGC 1851, where we see cluster 
stars lying along an "anomalous", less populated RGB (see e.g. 
Fig- SI- In this cluster there have been (controversial) claims of 
variations of the sum of C+N+0; this will be discussed in a pa- 
per in preparation. Here we simply propose that these stars could 
be the ones which are both C and N-rich, while the ones which 
have lower N abundance fall in with the bulk of the stars, even if 
C-rich. 

Finally, we considered another possible composition differ- 
ence; at the RGB bump t here is a change in the surface abun- 
dance of C and N (see e.g.. lGratton et alil2000h . i.e., C deer eases 
and N increases, due to a mixing episode. In C-normal, N-poor 
stars the increase in N after the RGB bump almost compensates 
the decrease in C, and the CN bands do not change their strength. 
Instead, in C-normal, N-rich stars, the N abundance is already 
high and the N increase has proportionally a smaller impact and 
does not compensate enough the decrease in C. Hence, there is 
a smaller difference between a N-rich and N-poor star if we take 
into account the evolutionary changes (see the correspondingly 
smaller differences in u and v and all correlated indices). 



Armed with these results, we may then examine the predic- 
tion of our spectral synthesis on the observed indices. 

nil = {v - b) - (b - y) is the classical metallicity index of 
the Stromgren photometry. In the present context, we notice that 
b and y (hence b - y, the classical temperature indicator) are in- 
sensitive to variations in the abundances of the CNO elements. 
For mi, the only dependence is then contained in the v mag- 
nitudes. Since the CN violet band is very weak in metal-poor 
and/or warm stars, nii separates "N-poor" and "N-rich" only for 
cool, metal-rich stars. 

ci - {u-v)-{v-b) is the classical gravity index of Stromgren 
photometry, measuring the Balmer jump. In the present context, 
it essentially measures the difference between the strength of the 
NH and CN violet bands (the latter actually weighted twice). 
Therefore it separates well "N-poor" and "N-rich" at low metal- 
licities (where the CN violet band is negligible), but it changes 
the sign of its sensitivity for cool metal-rich stars, where the CN 
violet band is quite strong. The same holds for Cy {- ci - (b-y)), 
which has the same sensitivity on CNO abundances of ci. 

Our new index 64 - (u - v) - (b - y) keeps constant the sign 
of its sensitivity to CNO abundances over a range of tempera- 
ture/metallicity/gravity larger than for mi and ci or c,., although 
it also fails for the coolest and most metal-rich stars. The reason 
of this more uniform behavior is that it weighs only once the v 
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Table 3. Model sensitivity of Stromgren photometric indices. 
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Notes, (a) A computed as C-normal,N-poor minus C-normal,N-rich, but with a composition changed after the RGB bump for C,N. 
C-normal,N-poor stars have: [C/Fe]=[N/Fe]=0, [O/Fe]=[Mg/Fe]=[Ca/Fe]=+0.4, [Ti/Fe]=0.2; 

C-normal,N-rich stars have [C/Fe]=-0.2, [N/Fe]=+1.3, [O/Fe]=-0.1 (and the same abundances for Mg, Ca, and Ti); 
C-rich,N-poor stars have [C/Fe]=+0.2, [N/Fe]=0.0, O and other elements as in first case; 
C-rich,N-rich stars have [C/Fe]=-l-0.2, [N/Fe]=+1.3, [O/Fe]=-0.1 and other elements as in the first case. 
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Fig. 12. (a) Synthetic spectra for a N-poor (magenta), a N-rich 
(black), and a C-rich star at the RGB bump of a GC with 
[Fe/H]=-1.23 (see Sect. 4 for the various abundances used). 
Over-imposed are the transmission curves for the uvbyCa filters; 
some features of interest are indicated (see text), (b) The same, 
but for a star above the RGB bump. 



band, and it is then less sensitive to the strength of the CN vio- 
let bands. In practice, the sensitivity of 64 on CNO abundances 
is all contained in the (m - v) term. The addition of the {b - y) 
term is however useful to reduce its sensitivity to temperature 
and reddening. In particular, since 5 a is practically insensitive 
to reddening, it can be safely used also for clusters with strong 
differential reddening. 

5. 1. Comparison with Sbordone et al. (2011) 

This very recent paper approached the problem of the photomet- 
ric signatures of multiple populations in GCs from the theoret- 
ical point of view. The approach and purposes were different, 
since the authors wanted to present there full blown isochrones 
and take the effect of helium in consideration as well, while we 
concentrated on the RGB. They then computed self-consistent 
stellar models, taking into account the different chemical compo- 
sition of first and second generations stars both for the stellar in- 
teriors and the stellar atmospheres. They computed models with 
several composition (standard mixture, NaCNO-enhanced, He- 
enhanced) and translated the theoretical isochrones into the ob- 
servational planes for Johnson and Stromgren photometry, show- 
ing the effects on the main sequence, SGB, and RGB. 

Our paper and theirs, while addressing the same problem, 
differ in several points; (i) while we have worked on the obser- 
vational data of several GCs, they did not present a direct com- 
parison with observed cases, deferring this to subsequent papers; 
(ii) we have not computed new, self-consistent atmospheres and 
fluxes, but only computed the variation in the absorption spec- 
trum due to a change in composition, using a standard model. 
Their approach is formally better; however the effect is small, as 
shown also by their paper, since the molecular bands present in 
these stars affect only a small fraction of the total flux and the 
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various models differ only in the most externa l parts, near the 
stellar surface. In fact, iGustafsson et alJ (l2008h considered the 
impact of mild and heavy CN cycling on the structure of the at- 
mospheres and found that it is negligible for the temperatures of 
interest in our computations (see their Fig. 7); (iii) in our paper 
the changes in composition due to the first dredge-up and the 
extra-mixing after the RGB bump are taken into consideration. 
This has significant impact on the derived results, since C and 
N are involved in these changes; (iv) perhaps most importantly, 
they do all calculations only for one value of metallicity, while 
we explore the entire range of metallicity of Galactic GCs. This 
is an important factor, since we have seen that there is a strong 
dependence on metallicity of the sensitivity of the various colors 
and indices to different compositions (see Tableland Fig.[T3ll. 

Given all these differences, it is legitimate to ask whether we 
are obtaining similar results. Taking our case at [Fe/H]=-1.63 
to compare to their -1.62 and looking at their Fig. 8, we obtain 
similar, but smaller, differences between N-rich and poor com- 
positions in u (and of course b and y, where N has no influence). 
There are instead discrepancies in the v filter, where we find dif- 
ferences and they do not. However, there is a large difference in 
C abundance between the two cases: we use [C/Fe]=-0.2, com- 
pared to their -0.6, meaning that their model has [C/H]= -2.23, 
i.e., is more similar to our [Fe/H]= -2.23 case, where we too do 
not see any effect in the v filter (the CN bands are really very 
weak in these models). 

The large difference in C abundance and the single, quite low 
metallicity could also explain the fact that their models predict 
larger differences in c,, between the different compositions (see 
their Fig. 15) than our computations (Table O. In practice, we 
have stronger CN bands in our spectra and this has an effect on 
the V filter, that is actually weighted twice in the definition of 
Cv, and outweighs the effect in the u filter, where the NH band 
dominates. A more complete comparison between the two works 
is difficult and outside the goals of the present work. It could 
be interesting, however, to merge the two approaches, to better 
define the properties of multiple populations in GCs. 

6. Summary 

In this paper we explored the capabilities of the Stromgren pho- 
tometric system to characterize multiple populations that we 
firmly know (from spectroscopy) to be harbored in all GCs. 
Thanks in particular to the u filter, which encompasses the 
molecular NH band at 3400 A, and the v filter, at least for the 
metal-rich GCs, it is possible to well separate N-poor (Na-poor, 
O-rich, first generation) stars from N-rich (Na-rich, O-poor, sec- 
ond generation) stars using Stromgren multicolor photometry. 

In our systematic exploration of the parameter space, us- 
ing Grundahl's (and coworkers) photometry made available by 
Calamida et al. (2007) for 9 GCs, we found that the classical 
metallicity index m\ is able, in the present context, to separate in 
N-strength classes only cool and metal-rich stars. 

Both the typical gravity indicator ci and its improved version 
Cy defined by Yong et al. (2008), work well in separating N-rich 
and N-poor stars at low metallicity, but give ambiguous results 
in the metal-rich and cool stars regime. 

We defined a new indicator, S^, that seems to works better 
than the previous indices over a larger range of metallicity, al- 
though even in this case the separation of different populations 
is not achieved well in the cool and metal-rich regime. 

The variations that are sampled by all these indices are well 
accounted for by changes in the abundances of CNO elements, as 
demonstrated by a number of comparison with synthetic spectra 



convolved with the transmission filter. However, further work 
is necessary, combining for instance our observational approach 
with the theoretical one chosen by Sbordone et al. (2011). 

In the present paper we have limited our analysis to the RGB 
stars to ensure the best photometric data and to compare directly 
with the available firm constraints given by the chemical abun- 
dances (e.g., the Na-O anticorrelation). The next step will be 
the extensive analysis of photometric data over wide fields and 
reaching also extremely blue (faint) HB and main sequence stars 
with comparable precision. The gain in statistics will permit to 
better study hot issues such as e.g., the radial distribution of mul- 
tiple stellar populations in GCs, which is tightly connected with 
the formation and evolution of these stellar aggregates. 
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